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ABSTRACT 

A  comprehensive  computer  program  for  designing  and 
evaluating  multidisciplinaiy  aircraft  flight  control  systems  is  de¬ 
scribed.  The  FlighrCAD  program  contains  a  variety  of  model¬ 
ing.  synthesis,  analysis,  and  simulation  alternatives.  The  pro¬ 
gram  will  be  used  to  produce  a  control  design  for  the  t99l 
AJAA  Controls  Design  Challenge  (to  be  completed  in  1992). 
FlightCAD  ultimately  will  implement  a  number  of  control  design 
techniques;  here  it  uses  direct  digital  synthesis  to  produce  a  ro¬ 
bust,  propordonal-integral-filter  controller  \kith  scheduled  linear- 
quadratic-Gaussian  gains  and  command  generator  tracking  of 
pilot  inputs.  The  FlightCAD  design  approach  is  reviewed,  and 
a  status  repon  is  presented. 

INTRODUCTION 

Designing  flight  control  systems  for  modem  aircraft  contin¬ 
ues  to  be  a  challenging  task;  high-performance  aircraft  are  be¬ 
coming  harder  to  control  as  a  consequence  of  lighter  weight, 
broader  flight  envelopes,  and  increased  performance  require¬ 
ments.  There  has  been  a  transition  from  inherently  stable  con¬ 
figurations  with  relatively  rigid  airframes  and  traditional  control 
surfaces  to  flexible,  unstable  configurations  with  multi-function, 
redundant  control  surfaces.  Flight  control  systems  with  me¬ 
chanical  linkages  and  limited-authority,  analog,  stability  aug¬ 
mentation  have  been  replaced  by  full-authority,  fly-by-wire, 
digital,  command  augmentation.  The  desire  for  "Level  1"  flying 
qualities  throughout  the  flight  envelope  leads  to  coupled,  high- 
gain  controllers  that  excite  structural  modes  and  interact  with 
control-actuator  dynamics,  and  there  is  increasing  need  to  inte¬ 
grate  flight  controls  with  engine  controls  and  load-alleviation 
functions.  These  are  imponant  precursors  to  multidisciplinary 
design,  in  which  aerodynamic,  structural,  propulsive,  and  con¬ 
trol  functions  are  considered  together. 

Currently  available  computer-aided  control  system  design 
(CACSD)  programs  facilitate  design  in  many  ways,  but  they  fall 
shon  in  flight-control-specific  features.  Factors  that  must  ^  ad¬ 
dressed  include; 

•  flight  dynamic  modeling 

•  approximating  aerodynamic,  structural, 
propulsive,  and  control  functions 

•  generating  nominal  flight  paths  and  flight  envelope 

•  satisfying  flying  qualities  criteria 

•  adapting  to  flight  condition 

•  analyzing  stability  and  performance  robustness 

These  capabilities  should  be  accessible  in  an  easily  understood, 
interactive  format  that  is  intuitive,  comprehen.sive,  and  requires 
little  supplemental  programming. 

The  FlightCAD  computer  program  is  intended  to  address 
these  needs.  It  will  have  broad  capabilities  for  modeling  aircraft 
systems  and  subsystems,  for  integrating  their  coupled  effects  in 
simulation,  and  for  designing  flight  control  systems  over  the  en¬ 
tire  flight  envelope.  The  program  uses  a  desktop  metaphor  to 
organize  the  design  process,  with  a  menu  bar,  pull-dow  n  lists  of 
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design  functions  and  alternatives,  dialog  boxes,  and  multiple 
display  windows.  FlightCAD  is  implemented  using  features  of 
the  NeXT  Computer  for  designing  user  interfaces,  integrating 
code  produced  in  several  programming  languages,  and  multi¬ 
tasking  within  a  UNIX  environment. 

In  the  remainder  of  the  paper,  the  1991  AlAA  Controls 
Design  Challenge  is  briefly  reviewed,  the  architecture  of 
FlightCAD  is  outlined,  the  Challenge  Design  control  logic  is 
presented,  and  the  current  status  of  development  is  outlined. 

1991  AlAA  CONTROLS  DESIGN  CHALLENGE 

The  1991  AlAA  Controls  Design  Challenge  presents  six- 
degree -of-freedom  dynamic,  aerodynamic,  and  thrust  data  for  a 
high-performance  aircraft  in  the  form  of  FORTRAN  code  and 
'vritten  specifications  |1].  Its  specified  goals  are  to  design  an 
automatic  digital  controller  that  a)  can  maintain  straight-and-level 
flight  at  a  specified  altitude  and  Mach  number,  b)  can  control 
normal  acceleration,  altitude,  and  Mach  number  in  a  constant-g 
turn,  and  c)  can  provide  a  level  acceleration  from  subsonic  to 
supersonic  flight.  The  Challenge  is  motivated  by  many  of  the 
factors  mentioned  above,  as  well  as  by  practical  issues  often 
overlooked  in  theoretical  development.  Control  design  solutions 
must  account  for  significant  nonlinearities,  time  variations,  and 
uncertainties  in  aircraft  dynamics;  limited  and  imprecise  sensors; 
and  control  surface  dynamics,  including  displacement  and  rate 
limits.  Bias  errors,  scale  factor  variations,  disturbances,  and 
noise  effects  also  must  be  considered. 

OVERVIEW  OF  FlightCAD 

FlightCAD  is  a  CACSD  program  focused  on  aircraft  dy¬ 
namics,  flight  control  systems,  stability,  and  performance.  It  is 
built  around  a  desktop  metaphor  that  features  pull-down  menus 
and  dialog  boxes  containing  design  alternatives,  and  documents 
with  multiple  layers  of  information.  The  program  runs  in  the 
NeXTStep  2.0  environment  on  a  NeXTCube  68040  Computer. 

On  first  stanup,  FlightCAD  presents  two  blank  windows, 
icons  of  frequently  used  software  "tools,"  and  a  menu  side-bar 
(Fig.  1).  One  window  contains  a  FlightCAD  document  (a 
CADdoc),  and  the  other  suppons  a  simple  text  editor  for  func¬ 
tion  entry  and  revision.  The  standard  NeXT  icon  dock  is  modi¬ 
fied  to  display  only  programs  relevant  to  FlightCAD.  The  top- 
level  menu  items  describe  FlightCAD's  functions  and  features. 
Several  of  the  items  are  common  to  most  window-based  appli¬ 
cations  (such  as  Info  and  Edith  those  that  are  specific  to 
FlightCAD  include; 

Block  Diagram  -  Define  and  display  the  structure  and  sub¬ 
structures  of  the  system's  elements,  as  depicted  in  the 
CADdoc. 

Model  -  Build  differential  equations,  algebraic  equations,  and 
transfer  functions  that  describe  the  reference  frame,  aircraft 
(including  inertial,  aerodynamic,  propulsive,  elastic 
properties),  actuators,  sensors,  and  flight  computer. 

Design  -  Compute  operating  points,  trajectories,  and  flight 
envelope.  Portray  nonlinear  functions.  Synthesize  control 
and  estimation  logic. 

Analv7»‘  -  Compute  time  end  frequency  responses.  Evaluate 
stability  and  performance.  Conduct  statistical  studies  of 
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Figure  1.  FlightCAD  at  first  startup. 


computational  results,  including  Stochastic  Robustness 
Analysis. 

Activate  -  List  of  all  open  documents.  Selecting  a  name  brings 
that  document  into  the  foreground. 

Menus  and  submenus  are  arranged  in  a  natural  hierarchy. 
In  the  example  (Fig.  2),  "Model"  is  the  top-level  selection, 
"Dynamics"  is  the  second-level  selection,  "Reference  Frame"  is 
the  third-level  selection,  and  "Flat  Eanh"  has  been  selected  as  the 
dynamic  modeling  frame  of  reference.  This  "point-and-click" 
process  defines  the  equations  of  motion  to  be  used  for  control 
design  and  simulation.  Most  remaining  choices  can  be  selected 
in  the  same  way,  and  the  NeXTStep  environment  makes  adding 
menu  items  an  easy  process. 

FlighiCAD's  focal  point  is  the  block  diagram.  The 
FlightCAD  document  contains  diagrams  drawn  by  the  user  to 
represent  the  system.  Each  block  has  several  labels  that  describe 
its  properties  (Fig.  3);  for  example,  a  given  block  might  have  a 
label,  its  contents  might  be  stored  in  a  filename,  its  function  may 
assume  continuous-  or  discrete-time  modeling  (‘CT’  or  *DT’), 
and  it  might  have  a  number  of  inputs  and  outputs.  Unknown 
properties  are  designated  with  a  question  mark.  A  block  can 
contain  a  function  (and  be  marked  with  an  ‘F’).  which  can  be 
edited  in  a  text-editor  window  (Fig.  4),  or  another  block  diapam 
(and  be  marked  with  a  ‘B’),  which  can  be  opened  and  examined 
in  a  separate  CADdoc.  In  this  way,  the  user  can  navigate 
through  the  existing  hierarchy  of  inner  and  outer  loops.  Blocks 
can  be  moved,  examined  and  connected  with  multi-variable  links 
using  the  menu  items.  Connections  between  blocks  can  carry  as 
many  variables  as  the  blcckisi  have  inputs  and  outputs.  All 
block  parameters  are  easily  changed  through  the  appropriate 
menu  items.  The  CADdoc  block  diagram  can  be  exponed  as  a 
PostSenpt  picture  for  inclusion  in  a  written  document. 


Block  diagram  templates  facilitate  system  development. 
Figure  3  represents  a  typical  flight  control  problem;  it  contains 
blocks  for  actuators,  sensors,  aircraft  dynamics,  the  environ¬ 
ment,  and  the  control  logic.  This  format  allows  quick  develop¬ 
ment  in  as  much  detail  as  desired.  If  a  given  block  is  not  re¬ 
quired.  its  function  is  set  to  unity;  hence,  a  control  designer 
could  simulate  just  the  plant  dynamics  at  first,  then  check  the  ef¬ 
fect  of  pressure  or  temperature  variations,  add  actuator  and  sen¬ 
sor  dynamics,  and  finish  with  a  controller  that  takes  all  of  this 
into  account. 

FlightCAD  translates  the  functions  represented  by  block  di¬ 
agrams  into  C  source  code.  At  each  step.  FlightCAD  incremen¬ 
tally  compiles  only  those  portions  of  the  code  that  have  changed 
since  the  last  compilation,  allowing  development  to  occur  in 
stages.  This  feature  supports  rapid  prototyping  and  comparison 
of  competing  designs. 

FlightCAD  contains  advanced  iteration  and  search  capabili¬ 
ties  that  enhance  modeling,  design,  and  analysis.  In  any  phase, 
the  controller  can  be  evaluated  at  selected,  tabulated,  or  random 
points  of  the  operating  range  space,  supponing  point  designs, 
designs  along  (or  in  the  vicinity  of)  nominal  path  histories,  or 
designs  that  span  the  entire  flight  envelope.  This  feature  is  use¬ 
ful  not  only  during  design  but  in  the  evaluation  phase,  when  the 
likelihood  of  satisfactory  stability  and  performance  must  be  de¬ 
termined. 

Control  and  estimation  design  algorithms  ultimately  will  in¬ 
clude  a  wide  range  of  alternatives,  from  classical  methods  for 
single-inpui/single-output  systems  to  modem  methods  of  multi- 
input/multi-output  design;  from  linear,  time-invariant  models  to 
nonlinear,  time-varying  models;  and  for  continuous  and  sam- 
pled-data  controllers.  (5ur  response  to  the  1991  AlAA  Controls 
Design  Challenge  focuses  on  gain-scheduled,  linear-quadratic- 
Gaussian  theory;  hence,  the  associated  elements  of  FlightCAD 
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Figure  2.  FlightCAD's  Menu  Smicture. 

will  be  developed  first.  Succeeding  efforts  will  broaden  the 
program  s  capabilities. 
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Our  Challenge  Design  has  a  conventional  block-diagram 
structure,  as  indicated  by  Figure  3.  The  Controller  block  con¬ 
tains  forward-  and  feedback-loop  elements,  and  the  Estimator 
block  contains  filters  for  model-based  noise  reduction. 
Feedback  elements  account  for  the  effects  of  low-frequency 
disturbances,  high-frequency  noise,  and  plant  uncertainty  using 
Proportional-Integral-Filter  (PIF)  compensation.  Forward-loop 
elements  generate  the  desired  state  trajectory  and  corresponding 
nominal  control  settings  using  a  Command-Generator-Tracker 
(CGT)  structure.  The  digital  flight  control  system  (DFCS)  takes 
an  incremental  form  to  facilitate  initialization,  interfacing  with  the 
nonlinear  environment,  and  switching  between  command 
modes  The  Estimator  block  contains  a  bank  of  reduced-order 
filters,  with  low-pass,  complementary,  and  notch  characteristics. 


Figure  4.  FlightCAD’s  Editor  Windows. 

as  appropriate.  These  filters  follow  a  Kalman  format,  although 
the  implementations  are  hierarchical  and  sub-optimal. 

The  Challenge  Design  DFCS  is  based  on  an  evolution  that 
began  with  the  study  of  coupled,  linear,  time-invariant  dynamic 
models  of  a  rapidly  maneuvering  aircraft  [2,3]  and  continued 
through  the  development  of  numerous  linear-quadraiic-Gaussian 
(LQG)  flight  control  designs  [4-19],  seven  of  which  were  tested 
in  flight.  The  dynamics  and  control  of  high-performance  aircraft 
were  investigated  in  |2-7].  culminating  in  suggested  designs  for 
Type  0  and  1  controllers  for  an  operational  Naval  aircraft.  The 
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Figure  3.  FlightCAD’s  Block  Diagram  Template. 


PIF  control  structure  was  developed  for  NASA's  VTOL 
Approach  and  Landing  Technology  (VALT  )  Program  [8-101, 
resulting  in  a  successful  flight  test  of  adaptive  attitude-  and  ve¬ 
locity-vector  command/control  laws  on  an  experimental  CH-47B 
helicopter.  Six  digital  control  laws  were  flight-tested  on 
Princeton  University's  two  Navion  variable-stability  aircraft, 
demonstrating  proportional,  proportional-filter  (PF),  and  pro¬ 
portional-integral  (PI)  linear-quadratic  control  [11,12],  direct 
side-force  control  via  implicit  model  following  [13],  adaptive 
lateral-directional  control  to  eliminate  "wing  rock"  in  fully  stalled 
flight  [14],  and  PIF  command  augmentation  [15,16].  PIF  con¬ 
trollers  were  designed  for  the  B-737  twin-jet  transport  used  in 
NASA’s  Advanced  Transport  Operating  Systems  (ATOPS) 
Program  [17,18].  The  PIF/CGT  approach  also  has  been  applied 
in  a  design  study  for  the  X-29  Forward  Swept-Wing 
Demonstrator  Aircr^t  [19]. 

The  Challenge  Design  extends  prior  developments  in  sev¬ 
eral  ways.  It  uses  the  PIF/CGT/LQG  model  to  produce  fully 
coupled  control  of  the  Challenge  Aircraft  throughout  the  design 
envelope.  Control  and  estimation  gains  at  each  of  several  fixed 
operating  points  will  be  chosen  to  satisfy  military  flying  qualities 
specifications  with  considerable  margins  in  the  initial  design. 
These  gains  then  will  be  adjusted  for  maximum  stability  and  per¬ 
formance  robustness  using  Stochastic  Robustness  Analysis  [20- 
22].  Artificial  neural  networks  will  be  used  to  adapt  control  and 
estimation  gains  to  changing  flight  condition  [23].  The 
Challenge  Design  will  be  developed  in  the  FlightCAD  environ¬ 
ment. 

Dynamic  Modeling 

The  aircraft  equations  of  motion  are  nonlinear  ordinary  dif¬ 
ferential  equations, 

x(t)  =  flx(0.u{f),w(t),pl  (1 

where  x(f)  is  the  n-dimensional  state,  u(t)  is  the  m-dimensional 
control,  w(f)  is  an  s-dimensional  disturbance,  and  p  is  an  /-vec¬ 
tor  of  parameters.  The  rigid-body  state  for  this  problem  is  de¬ 
scribed  by 

\~[\'  a  q  6  h  0  r  p  (p  y/]T  (2 

where  (V',  a.  P)  are  the  airspeed,  angle  of  attack,  and  sideslip 
angle  components  of  the  air-relative  velocity  vector,  {p.  q,  r)  are 
components  of  the  body-axis  angular-rate  vector,  ((p.  6.  if)  are 
the  Euler  angles  describing  roll,  pitch,  and  yaw  attitude,  and  h  is 
the  altitude.  The  first  five  components  of  x  are  called  longitu¬ 
dinal  variables,  and  the  second  five  represent  lateral-direc¬ 
tional  variables.  The  available  control  effectors  are  contained 
in. 

u  =  [5£  ST  SA  SR]^  (3 

corresponding  to  symmetric  deflection  of  the  left  and  right 
elevons  (or  stabilators),  thrust  command  (to  two  engines), 
asymmetric  deflection  of  left  and  right  ailerons,  and  rudder  de¬ 
flection.  The  first  two  control  elements  have  greatest  effect  on 
longitudinal  variables,  while  the  second  two  are  principally  lat¬ 
eral-directional  controls.  Wind  provides  the  principal  distur¬ 
bance  w;  it  can  be  expressed  as  perturbations  to  (V^  a.  p). 

The  state  (as  well  as  control  and  disturbance)  is  observed 
through  the  measurement  r-vector, 

z(t)  =  yU)  +  n(t)  =  h[x(r).u(r),w(r)]  -i-  n(r)  (4 

where  n(t)  is  an  r-dimensional  measurement-error  vector.  The 
measurement-vector  components  are, 

y  =  [as  at  as  p  q  r  \ff  e  0  hi  M  pT  Pa  a]  (5 


(as,  ai,  as)  are  orthogonal  components  of  acceleration  (not 
measured  at  the  center  of  mass),  (y/,  6.  0)  are  Euler  attitude  an¬ 
gles,  hi  is  inenial  altitude,  M  is  Mach  number,  and  (pr.  Pa) 
total  and  ambient  pressures. 

Linearized  models  find  widespread  use  in  control  design. 
Along  a  nominal  trajectory  specified  by  Xo(r),  Uo(t),  Wo(r),  and 
no(/)  for  t  in  (to,  tf),  perturbations  of  the  state  and  observation 
vectors  are  governed  approximately  by  linear,  time-varying 
equations, 

Ax(/)  =  F{t)Ax(t)  -f  G(t)Au(r)  -i-  L(r)Aw(f)  (6 

Az(/)  =  H„(t)Ax(r)  -*•  Hu(t)Au(r)  +  Hw(t)Aw(t)  -^  n(r)  (7 

F,  G,  L,  Hx,  Hu,  and  Hw  are  conformable  Jacobian  matrices 
expressing  sensitivities  to  the  penurbation  variables,  evaluated 
along  the  nominal  trajectory. 

F,  G,  and  L  panition  into  direct  and  coupling  blocks.  The 
upper-left  blocks  express  longitudinal  effects  on  longitudinal 
variables,  the  lower-right  blocks  express  lateral-directional  ef¬ 
fects  on  lateral-directional  variables,  and  the  off-diagonal  blocks 
induce  coupling  between  the  two  sets.  When  the  wing  span  is 
horizontal,  the  coupling  blocks  of  an  aerodynamically  symmetric 
aircraft  are  small,  as  non-zero  terms  arise  only  from  engine  gy¬ 
roscopic  effects.  In  turning  or  rolling  flight,  longitudinal  and 
lateral-directional  motions  are  coupled,  so  the  coupling  blocks  of 
F,  G,  and  L  may  have  significant  effect  In  a  steady  pullup,  the 
motions  are  not  coupled,  but  the  lateral-directional  blocks  are  al¬ 
tered  from  their  cruising-flight  values 

The  state  can  be  propagated  between  discrete  instants  of 
time  (fk.  fk+i,  ...  )  by, 

Axk+i  =  «l>kAxk  +  FkAuk  AkAwk  (8 

and  the  corresponding  measurement  penurbaiions  are 

Azk  =  Hx,,Axk  +  Hu^Auk  +  Hw^Awk  -h  nk  (9 

The  subscript  "k"  indicates  evaluation  at  fk-  Here,  O,  F,  and  A 
have  the  same  dimensions  as  F,  G,  and  L  and  are  derived  from 
the  system's  state  transition  propenies  for  a  given  sampling  in¬ 
terval,  A/  =  fk+i  -  fk-  For  present  purposes,  a  constant  control 
sampling  rate  (=  1/Af)  of  50  per  sec  appears  appropriate,  while 
gain  schedules  may  be  updated  at  a  slower  rate. 

Dynamic  Compensation 

Control  and  estimation  logic  for  the  nonlinear  plant  (eq.  1 
and  2)  can  be  expressed  as  a  dynamic  compensator. 


Auk  =  -Ck^k  (10 

^k+l  =  ^k4k  +  ©kAuk  +  Kk[zk  -  h(Xk,Uk)]  (11 

^k  =  lAxkT  Xk'^F  (12 

Uk  =  Uo,  -t-Auk  (13 

k 

Xk  =  Xo.  +Axk  (14 

k 


This  linear,  time-varying  structure  is  equivalent  to  a  feedback 
control  law  (eq.  10)  operating  on  the  internal  state  estimate  Ax 
contained  in  the  (n  -t-  )t)-dimensional  4k  (eq-  12)-  Xk  is  a  Jt-vec- 
tor  of  compensation  components  that  may  include  integrals  of 
state  elements  or  an  explicit  command  model  to  be  followed.  Ck 
and  Kk  are  selected  to  provide  satisfactory  nominal  response, 
and  they  may  vary  in  time.  Vk  and  ©k  include  nominal  values 
of  Ok  and  Fk  plus  integrating  (i.e.,  accumulating)  or  filtering 
operations  associated  with  Xk  Ck,  Kk.  'Fk,  and  ©k  are  repre- 


samcd  as  full-order  matrices  here,  but  they  may  be  panitioned, 
reduced,  and  simplified  in  DFCS  implementation. 

The  desired  state  and  corresponding  control  for  the  nonlin¬ 
ear  plant,  Xo^and  Uoj^,  enter  as  in  eq.  13  and  14,  and  they  can  be 

produced  by  the  Command  Generator.  The  nominal  control  set¬ 
ting  is  quite  sensitive  to  parameter  varianons  and  disturbance  in¬ 
put.  The  incremental  form  control  law  eliminates  the  need  to 
know  Uo^  by  computing  the  present  control  command  as  a  per¬ 
turbation  to  the  previous  command,  which  implicitly  contains  the 
nominal  control  value: 


Uk  =  Uk-l  -  (Cic  -  Ck-l  4k  i)  (15 

Proporlional-Integral-Filter/Command-Generalor 
Tracker  Control 

The  Proponional-lntegral-Fiher  (PIF)  structure  augments 
aircraft  dynamics  with  both  the  control  rate  dynamics  and  com¬ 
mand  output  integration.  It  allows  control  displacements  and 
rates  to  be  limited  implicitly,  and  it  provides  explicit  integration 
of  the  error  between  the  desired  and  actual  command  outputs. 
The  control  law  is  derived  using  the  line.arized  aircraft  model, 
augmented  with  equations  for  conffol  rate  and  command  output 
integration  [16]: 
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The  corresponding  sampled-data  form  is. 
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Equation  16  requires  that  Ax  and  Au  be  integrated  to  form 
A^  in  eq.  17,  and  a  first-order  data  hold  would  be  needed  to  ob¬ 
tain  AU|j  from  Av.  Euler  integration  is  used  in  both  cases,  and 

the  couplings  to  Ax  and  A^  are  eliminated,  yielding  the  follow¬ 
ing  model: 
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where  Ax  is  the  perturbation  from  the  optimal  state  history,  Au 

is  the  perturbation  from  the  optimal  control  history,  and  Av  is 
the  penurbation  from  the  optimal  control  rate  history. 

The  cost  function  to  be  minimized  by  control  is, 


The  desired  state  trajectory  can  be  propagated  by  a  com¬ 
mand  generator,  which  shapes  pilot  inputs  to  produce  a  desirable 
nominal  response.  In  combination  with  suitable  logic  for  regu¬ 
lation  (e  g.,  PiF),  the  controller  becomes  a  command-generator- 
trackcr  (CGT).  CGT  is  a  form  of  explicit  model  following, 
w  here  a  desired  set  of  dynamics  is  simulated  and  the  actual  sys¬ 
tem  is  controlled  as  a  function  of  the  error  between  the  two  [27], 
This  approach  allows  for  uncenainty  in  the  plant  dynamics 
(reflected  in  F  and  G)  and  typically  requires  high  controller 
gains  for  close  tracking  of  the  generated  commands.  However, 
a  combination  of  explicit  and  implicit  model  following  has  been 
shown  to  produce  satisfactory  controllers  with  modest  gains 
125).  For  the  initial  design,  command-generator  logic  will  be 
approximated  by  simpler  matrix  relationships  that  achieve  satis¬ 
factory  steady-state  response  [26],  Additional  details  of  the 
PIF/CGT  controller  can  be  found  in  [16,18], 

State  Estimation 

The  PIF/CGT  controller  is  designed  with  the  assumption  of 
full  state  feedback.  Because  all  state  components  are  not  ob¬ 
served  directly  and  all  measurements  are  subject  to  noise,  bias, 
and  scale  factor  error,  state  estimation  is  needed.  In  addition  to 
reducing  noise  effects  in  all  observations,  the  main  goals  are  to 
estimate  the  air-relative  velocity  components,  altitude,  and  atti¬ 
tude;  to  estimate  biases  and  scale  factors  as  required;  and  to 
provide  quantities  needed  for  gain  scheduling  (see  below). 

Prior  experience  [8-11]  suggests  that  several  reduceid-order 
filters  are  preferable  to  a  single  full-state  estimator.  The  prelimi¬ 
nary  estimator  suite  is  as  follows:  a)  Three  uncoupled  first-order 
angular-rate  filters,  b)  one  velocity  vector  filter  (with 
accelerometer  bias  estimate),  c)  one  altitude  filter  (with  altimeter 
bias  estimate),  d)  one  attitude  vector  filter  (possibly  using 
quaternions  internally),  and  e)  denvation  of  gain  scheduling 
parameters. 

Gain  Scheduling  via  Neural  Networks 

Challenge  Design  DFCS  gains  will  be  calculated  at  operat¬ 
ing  points  throughout  the  flight  envelope,  as  specified  by  angle 
of  attack,  Mach  number,  dynamic  pressure,  and  wind-axis  roll 
rate  (e.g.,  [2]).  Anificial  neural  networks  provide  an  exciting 
new  alternative  for  scheduling  these  gains  because  they  can  gen¬ 
eralize  nonlinear  functions  of  many  variables  with  little  or  no 
prior  knowledge  of  those  functions’  shapes.  Many  neural  net¬ 
works  receiving  current  attention  are  memory-less  expressions 
that  approximate  functions  of  the  form. 

(21 

where  x  and  y  are  input  and  output  vectors,  and  /(•)  is  the  rela¬ 
tionship  between  them.  Neural  networks  can  be  considered 
generalized  spline  functions  that  identify  efficient  input-output 
mappings  from  observations  [23].  Rather  than  approximating 
eq.  28  by  a  series,  table  look-up,  or  traditional  splines,  an  N- 
layer  neural  network  represents  the  function  by  recursive  opera¬ 
tions. 


where  the  weighting  matrices  are  chosen  to  meet  implicit  model¬ 
following  goals  associated  with  good  handling  qualities  [24,25], 
The  resulting  sampled-data  optimal  control  law  takes  the  form 


x(k)  =  s(k)[W(k-l)x0‘-')]  =  sO')[TiO')] ,  k  =  1  to  N  (22 

where  y  =  x(^)  and  x  =  xW.  W(*‘'')  is  a  matrix  of  weighting 
factors  determined  by  the  learning  process,  and  s('')[.]  is  an  acti¬ 
vation-function  vector  whose  elements  normally  are  identical, 
scalar,  nonlinear  functions  Oi(Tii)  appearing  at  each  network 
node: 

sOOjnO')]  =  [aidliC'))  ...On(tln<*>)]’^  (23 


Av^  ,  =  C|Ax  CjAu -t- C^A?  (20 

C|.  Cs.  and  Ct  are  the  solutions  to  the  algebraic  Riccati  equa¬ 
tions  associated  with  the  regulator  [26], 


One  of  the  inputs  to  each  layer  is  a  unity  threshold  element  that 
biases  the  activation-function  output. 

The  sigmoid  is  commonly  used  as  the  anificial  neuron, 
though  other  functions  such  as  the  derivative  of  the  sigmoid  or 
the  radial  basis  function  can  be  used.  The  sigmoid  is  a  saturat- 


ii;g  function  written  as  o(Ti)  =1/(1  +  e’’l)  for  output  in  (0,1). 
Learning  algorithms  for  the  elements  of  typically  involve 
either  a  gradient  search  or  an  extended  Kalman  filter. 

In  the  present  application,  each  gain  element  will  be  repre¬ 
sented  by  a  neural  network  with  four  inputs  (the  parameters 
specifying  the  operating  space  for  gain  calculations)  and  a  single 
hidden  layer  of  four  nodes.  The  networks  need  not  be  trained 
on-line;  hence,  they  do  not  add  an  undue  amount  of  complexiiy. 
They  do,  however,  provide  substantial  power  for  portraying 
gain  schedules  efficiently  and  with  small  error. 

Stochastic  Robustness  Analysis 

Robustness  of  the  Challenge  Design  DFCS  can  be  assessed 
in  a  number  of  ways.  Gain  and  phase  margins  provide  weak  in¬ 
formation  about  the  robustness  of  single-input/single-output 
systems.  They  are  not  reliable  indicators  of  the  likelihood  that 
stability  and  performance  goals  will  be  met  under  parameter 
uncertainty  (28]  because  parameter  variations  change  the  shape 
of  the  Nyquist  contour  as  well  as  its  gain  and  phase  [26|. 
Similar  notions  have  been  applied  to  the  mulri-input/multi-output 
case  through  singular-value  analysis,  but  the  same  Nyquist 
contour  effects  occur;  hence,  a  new  approach  to  defining  control 
system  robustness  is  needed. 

Stochastic  Robustness  Analysis  is  such  an  approach  [20- 
22|.  It  is  based  on  the  premise  that  each  parameter  variation  can 
be  modeled  by  a  probability  density  function,  which  may  be 
bounded,  non-Gaussian,  and  correlated  with  other  parameters. 
This  probability  description  can  be  used  to  shape  the  output  of  a 
random  number  generator  which,  taken  with  the  outputs  for  all 
other  uncertain  parameters,  forms  a  set  of  parameters  for  a  single 
evaluation  of  stability  and  performance  metrics.  The  metrics  are 
noted,  and  the  process  is  repeated  in  a  "Monte  Carlo"  analysis. 
As  the  number  of  trials  grows,  the  tested  parameter  distributions 
approach  their  models,  and  the  probability  distributions  of  the 
stability  and  performance  metrics  are  estimated  within  com¬ 
putable  confidence  intervals.  These  estimates  are  compared  with 
target  values  to  determine  whether  or  not  the  closed-loop  system 
possesses  adequate  robustness. 

In  the  present  application,  we  are  particularly  concerned 
with  maintaining  a  vanishingly  small  probability  of  closed-loop 
instability,  a  high  probability  of  satisfactory  performance,  and  a 
reasonable  level  of  assurance  of  low  turbulence  response.  The 
appropriate  military  specifications  for  flying  qualities  and  flight 
control  systems  will  establish  target  levels  for  these  evaluations. 

STATUS  OF  THE  PROJECT 

The  AlAA  Design  Challenge  model  has  been  translated 
from  FORTRAN  to  C  code  and  has  been  compiled  on  the  NeXT 
Computer.  Results  from  the  C  version  are  being  compared  with 
those  generated  by  the  FORTRAN  version.  All  further  devel¬ 
opment  will  be  based  on  the  C  version. 

The  FlightCAD  program  is  operational  in  part;  the  user  in¬ 
terface  is  complete,  and  some  of  the  simulation  routines  and 
f,j..ctions  are  working.  The  block  diagram  manipulation  code  is 
not  complete  and  is  the  major  focus  of  current  work. 
FlightCAD  can  manipulate  text  documents  and  compile  them, 
and  can  open  completed  CADdocs.  Block  connection  and 
compilation  also  is  on  the  near-term  agenda. 

CONCLUSION 

A  flight  control  design  with  proportional-integral-filter 
structure,  command-generator-tracking,  linear-quadratic- 
Gaussian  gain  calculation,  and  gain  scheduling  is  proposed  as  a 
solution  to  the  AlAA  Controls  Design  Challenge.  The  design 
has  practical  engineering  sipificance,  making  it  suitable  for 
flight  test  and  future  application.  The  control  design  will  be  de¬ 
veloped  with  FlightCAD.  a  computer  program  that  presents  new 
paradigms  for  the  multidisciplinary  design  of  future  flight  con¬ 
trol  systems. 
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